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ABSTRACT 

We measured the radii of 7 low and very low-mass stars using long baseline interferometry with the VLTI interferometer and its VINCI 
and AMBER near-infrared recombiners. We use these new data, together with literature measurements, to examine the luminosity- 
radius and mass-radius relations for K and M dwarfs. The precision of the new interferometric radii now competes with what can be 
obtained for double-lined eclipsing binaries. Interferometry provides access to much less active stars, as well as to stars with much 
better measured distances and luminosities, and therefore complements the information obtained from eclipsing systems. The radii of 
magnetically quiet late-K to M dwarfs match the predictions of stellar evolution models very well, providing direct confirmation that 
magnetic activity explains the discrepancy that was recently found for magnetically active eclipsing systems. The radii of the early K 
dwarfs are well reproduced for a mixing length parameter that approaches the solar value, as qualitatively expected. 
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1. Introduction 

Measuring accurate masses, radii, and luminosities, for low 
and very low-mass stars has always been observationally chal- 
lenging. Precise individual masses and radii can be measured 
through combined photometric and spectroscopic observations 
of double-lined eclipsing binaries, but to date only ~15 such sys- 
tems are known with masses under 1 solar mass, and a few of 
those are too distant and faint for high precision work. Perhaps 
even more importantly, low mass eclipsing binaries tend, al- 
most by construction, to be fast rotators and hence magnetically 
veiy active. Evolutionary models have been found to system- 
atically underestimate the radii of very low mass eclipsing bi- 
naries ( [Torres & Ribas|[2002j l, and the uniformally high activ- 
ity level of these objects is currently the leading explanation for 
that discrepancy. 'Chabrier et al. ( |2007) l find that increased sur- 
face spots coverage, and for partly convective stars convection 
quenching by strong magnetic fields, can inflate the stellar radius 
by amounts which qualitatively match the observed discrepancy. 

Direct tests of that prediction, and validation of the 1-D 
structural models on objects which better match their assump- 
tions, need radius measurements for slowly rotating and mag- 
netically quiet very low-mass stars. Strong observational selec- 
tion effects unfortunately ensure that all known eclipsing bina- 
ries have sufficiently short orbital periods that they are tidally 
synchronised, and therefore in turn that all are fast rotators. 
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Measurements of magnetically quiet slow rotators therefore have 
to use a different observing technique, long-baseline optical or 
infrared interferometry. [Lane et al. (2001 1 and Segransan et al. 
( [2003| l both demonstrated 1-5% radius precision for stars that 
are only partially resolved. 

The mass of these single, isolated, stars is not directly ac- 
cessible, and, strictly speaking, both [Lane et "aT] ( [200 1) ) and 
Segransan et al. ( 2003| l therefore probe the luminosity-radius re- 
lation rather than the mass-radius one. Accurate mass and lumi- 
nosity measurements for M dwarfs (e.g. Segransan et al.[[2000 



however demonstrate that their K-band mass-luminosity relation 
has very low dispersion ( Delfosse et al. 20001), and mass is there- 



fore largely interchangeable with absolute K-band luminosity. 

Here we present direct angular diameter measurements 
for seven single K0.5 to M5.5 dwarfs, obtained with the 
VINCI and AMBER instruments on the Very Large Telescope 
Interferometer (VLTI) between 2003 and 2008. Section 2 de- 
scribes those observations, the data analysis, and the angular di- 
ameter determination. Section 3 discusses the luminosity-radius 
and mass-radius relation for very low-mass stars in the light of 
the new measurements, and compares the empirical relations 
with theoretical predictions. 

2. Observations and data analysis 

2.1. Sample 

The target list was largely determined by the capabilities of the 
VLTI at the time of the observations. The limiting correlated 
magnitude of the low spectral resolution mode of AMBER on 
the 1.8m auxiliary telescopes was K-4 in 2007, and improved 
to K=5.5 in 2008. We consequently selected targets with ap- 
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parent magnitudes between K=2.18 and K=4.38. Ongoing im- 
provements to the VLTI infrastucture are expected to provide 
access to fainter targets. We also computed the expected angular 
diameters from flux-colour relations and literature photometry, 
and selected targets with a predicted diameter above 0.9 mas. 
This translates into a visibility of at most 0.8 on a 128 m base- 
line in the H-band, as needed to measure diameters to a few % 
uncertainty with the current amplitude calibration precision of 
AMBER. One object, GJ 879, was observed as a backup target 
during an unrelated observing program, and does not fulfill this 
minimum angular diameter specification. 

2.2. Observations 

2.2.1. VINCI Observations 

GJ845 A (e Ind), GJ 166 A (DY Eri), GJ570A (KX Lib) and 
GJ663A were observed on the ESO Very Large Telescope 
Interferometer (VLTI) using its commissionning instrument, 
VINC I (IKervella et al.|2000| l with the two 35 cm test siderostats. 
Table [T summarises the observation details. VINCI operated 
in the K-band and used single-mode optical fiber couplers to 
recombine the light from two telescopes, and modulated the 
optical path difference around the white light fringe to produce 
interference fringes. This recombination scheme, first used 
in the FLUOR instrument ( |Coude Du Foresto et al.| 11998) , 
produces high precision visibility values, thanks to the efficient 
conditioning of the incoming wavefronts by the single mode 
fibers, to photometric monitoring of the light coupled into each 
input fibers, and to fast scanning of the high quality fringes. 
( Kervella et al.||2004a I extensively describe the data reduction 
ibr VINCI. 

Our observing strategy alternated sequences of several hun- 
dred fringe scans on the target star and on spatially close 
calibrator stars, to efliciently sample the temporal and spatial 
structure of the atmospheric and instrumental transfer function. 
Adherence to this strategy could unfortunately not always be 
strict, since scientific observations often had to give way to VLTI 
commissionning activities. For GJ 166 A, in particular, we could 
only keep three data points since all other measurements had no 
acceptably close calibrator observations. That target additionally 
was observed under poorer atmospheric conditions than all all 
other sources, and with two calibrators respectively located 78 
and 113 degrees away. It is by far our worst quality measure- 
ment, and is thus discarded from this study. 

2.2.2. AMBER Observations 

We used the AMBER (Astronomical Multi-BEam combineR) 
recombiner of the VLTI to measure the radiii of GJ166A, 
GJ887, Proximo (GJ551) and GJ879. Table [T] summarizes the 
observing circumstances. AMBER uses single-mode fibers for 
wavefront filtering and produces spectrally-dispersed fringes in 
the J, H, and K near-infrared bands ( |Petrov et al.|2007l l. 
We used the 1.8m-diameter VLT auxilliary telescopes (AT) on 
the AO-KO-Gl baseline triplet, which at the time of our obser- 
vations offered the longest available baselines and therefore the 
highest angular resolution. We selected the low spectral resolu- 
tion mode (low-JHK) of AMBER, which covers the J, H and 
K bands with R = 30, and adopted a 50 ms exposure time. 
AMBER unfortunately has poor J-band sensitivity, and we de- 
tected no fringes in that band. The H-band fringes, on the other 
hand, probe significantly higher spatial frequencies than the K- 



band VINCI would have on the same baselines. We mostly chose 
to not use the FINITO fringe tracker ( Le Bouquin et al.|2008| l, 
since our targets were at best close to the H=3 hmiting mag- 
nitude of FINITO. The fringes would thus not have been suf- 
ficiently stabilized to allow much longer integration times and 
compensate the 80% FINITO levy on the H-band flux. 
Since accurate absolute visibility calibration with AMBER had 
not been demonstrated when we planned the observations, we re- 
quested a very conservative observing strategy. Observations of 
up to four different amplitude calibrators were interleaved with 
those of each science target, to closely monitor the instrumen- 
tal and atmospheric transfer function. Those calibrators were 
chosen from Merand et al., (200 4) for angular proximity, well 
constrained predicted visibilities, and an approximate magnitude 
and color match to the science targets. The last requirement had 
to be relaxed somewhat for the M dwarf targets, since we could 
locate no apropriate M-type calibrators. For those stars we thus 
used K giant calibrators, which remain fairly close in near-IR 
colors. Since AMBER had never been used to measure precise 
angular radii, we chose to reobserve two stars previously mea- 
sured with VINCI, GJ 551 and GJ 887. VINCI has a well estab- 
lished record of accurate amplitude calibration, and those three 
stars provide a valuable check on any potential systematics in 
the AMBER measurements. 

GJ 879 was observed as a backup target in very poor atmo- 
spheric conditions (seeing FWHM > 1.5" and tq < 2ms), and 
has the smallest angular radius in our sample (highest V^). Its 
radius measurement, as a consequence, has significantly larger 
error bars than that of any other AMBER source. 

2.3. Data reduction 

With angular sizes under 2 mas, our targets are only partially 
resolved on the longest baselines available for our observations 
(128 m on the AO-KO basehne with AMBER, and 140m on the 
B3-M0 baseline with VINCI). Their squared visibilites re- 
main above 0.5 in the H-band. We therefore cannot derive their 
angular diameters from just the location of the first null of the 
visibility function, and instead need accurate calibration of the 
visibilities. The very partial resolution, on the other hand, en- 
sures that bandwidth-smearing effects are negligible. At our pre- 
cision level on the visibilities, accounting for the finite spectral 
bandwidth would only become necessary for y2 below 0.3 for 
VINCI, and well beyond the first visibility null for AMBER. 

The end products of both the VINCI and flie AMBER 
pipelines consist of coherence factors (ji^) together with an in- 
ternal error estimate on that quantity. The coherence factor is 
related to the squared visibility, V^, through: 



V 



(1) 



where is the squared transfer function of the instrument plus 
the atmosphere. 

Accurate calibration of the absolute squared visibilities 
therefore critically depends on a well understood transfer func- 
tion. That function is sensitive to both instrument stability and 
atmospheric conditions, and can fluctuate during a night. To as- 
sess its stability, we calculated the transfer function for every cal- 
ibrator exposure during a night, disregarding only those datasets 
for which too few scans/frames passed the reduction pipeline's 
thresholds to ensure the statistical significance of the resulting 



B.-O. Demory et al.: Mass-radius relation of low and very low-mass stars revisited with the VLTI 



3 



coherence factor measurement. In most cases, the transfer func- 
tion for each target measurement was evaluated from two dif- 
ferent calibrators, providing some control for temporarily dete- 
riorated atmospheric conditions. Using two calibrators also pro- 
tects against systematics introduced by a poorly chosen calibra- 
tor, such as unrecognised binaries. When no calibrator was ob- 
served immediately before or after a target point, we adopted 
the mean of the transfer function for the two nearest calibrators, 
provided they were observed within Ih of the science measure- 
ment. A few observations had to be discarded because no suffi- 
ciently close pair of calibrator observations was available. Table 
[2]and[3]respectively summarise the calibrators' properties for the 
AMBER and the VINCI observations. 

The contribution of the calibration to the visibility error bars 
accounts for both the statistical uncertainty on the coherence fac- 
tor and the propagation of the uncertainty on the calibrators' di- 
ameters. Since the different measurements for a given tar- 
get share any systematic error on the calibrators' diameters, they 
cannot be considered as fully independent. Those correlations 
are accounted for in the error bars, using the method described 
by Perrin ( 2003 ). That method is directly applicable to VINCI 
observations, but needed adaptation for AMBER, as we detail 
later 



2.3.1. VINCI data reduction 

We used V. 3.1 of the standard VINCI reduction pipeline 
( Kervella et al.| [200 4a). We used the wavelet spectral density 
( Segransan et al.|I999j l as a visibility estimator, as we found that 
it more robustly removes pistonned interferograms than the more 
common Fourier analysis. The commissionning state of the in- 
strument and of the VLTI array often affected the observing strat- 
egy, and we had to discard a significant number of measurements 
for which no calibrator observations where recorded at the same 
scan frequency. Kervella et al. ( 2004a| l showed, from observa- 
tions of brighter calibrators, that the squared transfer function 
of VINCI instrument is stable to within 1.5% during a night. We 
could therefore average the transfer function measurements from 
all calibrators observed during one night, and the observed dis- 
persion mainly reflects the statistical errors on these individual 
measurements. Those are the main source of uncertainties on our 
final radii obtained with VINCI. 



2.3.2. AMBER data reduction 

General description We used the AMBER data reduction 
pipeline described by Tatulli et al.| ( [2007[ ). After data reduction, 
we noted structures in the V'^ data obviously resulting from cor- 
relations in the dataset presumably due to non-optimal pipeline 
settings. Indeed, in february 2008, the AMBER TaskForce team 
( [Malbet et al. 2008 ) made a report of the displacement of the 
photometric channels with respect to the interferometric chan- 
nel between 2004 and 2008. Correlations occur on several pixels 
if those displacements are not correctly calibrated, a step that is 
achieved through the AMBER standard calibration matrix com- 
putation. During our October 2007 observing run, the spectro- 
graph entry slit was tilted and badly focused, thus propagating 
correlations over 4 to 5 pixels on the detector while AMBER 
spectral resolution is usually sampled over 2 pixels. We cor- 
rected this effect by taking into account integer pixel channel off- 
sets during spectral reshifting procedure, to avoid additional cor- 
relations to appear during the subtraction of the bad-pixel map 
at the sub-pixel level. None of our targets had sufficient J-band 



flux to offset the poor transmission of AMBER in that band. We 
therefore discarded the J-band data as well as wavelengths af- 
fected by telluric absorption, only keeping the centers of the H 
and K bands, 1.65 - 1.85 yum and 2.10 - 2.40 //m. 

We concatenated all observations of each object (usually 5 
sets of 1000 frames) into a single dataset from which we could 
more robustly select on fringe SNR, to reject pistonned interfero- 
grams as well as scans with low flux on one of the two telescopes 
in a pair. We verified that selecting the best 20% to 75% frames 
produced similar results, and therefore that the details of the se- 
lection do not unduly affect the outcome. Our final was to keep 
the best 20%, as a trade-off between fringe jitter suppression and 
increased noise. Selection on an absolute SNR threshold would 
more effectively reject blurred fringes (which bias down V^) than 
accepting a specified fraction of the data, but that alternate mode 
is not available in the current version of the AMBER pipeline. 
Another limitation of that data reduction package is that it uses a 
common threshold for all 3 baselines, in spite of their quite dif- 
ferent throughputs. Those limitations would become more criti- 
cal for datasets containing smaller number of scans than we used 
here. 



Transfer function As discussed above, a well understood trans- 
fer function is critical to measuring absolute visibilities. The sta- 
bility of the VINCI transfer function is well established, and a 
wealth of absolute visibilities have been published with that in- 
strument. The reliability of AMBER for absolute measure- 
ments, on the other hand, has been questioned and the stability 
of the AMBER needs to be established. Figure [l] shows the 
squared transfer functions of the 3 AMBER baselines during the 
night of 27 October 2007. The atmospheric conditions during the 
first half of the night were representative of Paranal, with a 0.8 
arcsec seeing and a coherence time slightly above 3 ms. With 
a 2 to 4% rms dispersion for both bands on baselines 1 and 3 
(90 and 128 m. length respectively) and in K-band for baseline 
2 (90 m.), AMBER approaches the stability of the VINCI trans- 
fer function. The end of the night had severely degraded atmo- 
spheric conditions (2.5 arcsec seeing and 1 ms coherence time), 
and the transfer function degraded only moderately, except on 
baseline 2, that had shown a 13% rms dispersion. 

Error bars computation The reliability of absolute mea- 
sured with AMBER has not yet been well established, and eval- 
uating realistic error bars therefore needs close attention. The 16 
to 18 spectral channels which we usually kept (7 to 8 in the H- 
band, and 9 to 10 in the K-band) are measured simultaneously. 
They therefore share the same atmosphere, as well as any er- 
ror on the angular diameter of the calibrator(s). If those factors 
dominate over statistical noise, the individual channels become 
highly non-independent. 

To quantify theses correlations between spectral channels, 
we generalise the formalism developed by Perrin (2003) for a 
two-channel combiner. For two Gaussian distributions of V- se- 
ries, the correlation coefficient between spectral channel k and a 
reference channel, r, is: 



((yj-yViyl-yl)) 



(2) 



we computed p^,„i( through Monte-Carlo simulations from the 
mean and variance of and V^^ . This provides a global corre- 
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Fig. 2. Calibrated squared visibilities from VINCI and best-fit 
LD disk model (solid) for GJ663A vs. spatial frequencies. 1-cr 
(dash) and 3-cr (dot) uncertainties are also indicated. 



Time [h] 

Fig. 1. AMBER squared transfer function for Oct 27, 2007. The 
upper, middle and bottom panel respectively represent baselines 
3, 2 and 1. Circles and squares respectively represent the values 
for the center wavelengths of the K-band (2.25/im) and of the H- 
band (1.76/im), The corresponding rms dispersions are reported 
in each panel. 



lation factor for a given band and to estimate amplitude of the 
error bar that is independent, and then to compute realistic er- 
ror on the final diameter. As expected, we found that data in 
a same spectral band are highly correlated and thus biases the 
final result if correlations are not taken into account. 

Error estimates on final radii obtained with AMBER mainly 
come from transfer function uncertainties and correlations be- 
tween squared visibilities. Those factors take more importance 
as the coherence time degrades. 

2.4. Data analysis 

2.4.1. Limb-darkened diameters 

At the level of accuracy achieved on the diameter determination 
for K dwarfs, discrepancy between uniform disk (UD) and limb- 
darkened disk (LD) is significant. Therefore, we used the non- 
linear limb-darkening law describing the intensity distribution of 
the star disk from |Claret| ( |2000| ) : 



k=l 



(3) 



where 1(1) is the specific intensity at the center of the disk, 
fi = cosy, J being the angle between the line of sight and 
the emergent intensity, and the limb-darkening coefficients. 
Tefj- and log g for each target are shown in table |4] with the 
corresponding references. We used those parameters, added to 
the photometric band (K for VINCI, H-hK for AMBER) and 
micro-turbulence velocity (assumed to be Yj - 2km/s), to select 



the corresponding limb-darkening coefficients for the PHOENIX 
models. 

Then, w e adju sted a limb-darkened disk that is given by 



Davis et al. (2000 1 to the data 



X' dpil (p)pih {nBOLD/A (l - IJ^-f'^) 
Jo dfil (fi) 



(4) 



We used the same LD coefficients for H and K bands since cor- 
responding discrepancy has been evaluated at the 0.05% level on 
the final radius, negligible as compared to the error bar amplitude 
obtained with AMBER. Uncertainties on the same coefficients 
due to a slightly different T^fj and log g can also be neglected, a 
200K change corresponding to a 0.01% on final radius estimate. 
Table|5]lists the derived UD and LD angular diameters and thek 
corresponding errors for both instruments. Figures [2] [3] and [AT| 
show data from VINCI with corresponding fitted LD mod- 
els for GJ663 A, GJ845 and GJ570A respectively. Figures [4] 



A.2 



and 



A.3 



show data from AMBER with correspond- 
ihgTitted LDlnodels for GJ 887, GJ 166 A, GJ 551 and GJ 879 
respectively. 



2.4.2. Instrument systematics 

Assessing systematics is essential to consider when reaching a 
few percent precision on angular diameters. We thus wanted to 
check for consistency between both instruments on GJ 887 and 
GJ 55 1 . The results we obtain are shown in Table |5] Angular 
diameter determinations are consistent for both instruments. 
GJ 887 has been observed in optimal conditions in both cases 
and the agreement is good at 1-cr level, as it is for GJ551. 
Bracketing each source point with two calibrators allow a better 
sampling of the transfer function, thus slightly reducing those ef- 
fects as shown on fig. [T] When good atmospheric conditions are 
met and a proper calibration applied, systematics on H and K 
band can be expected to be at 2% level, slightly above VINCI's. 
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Fig. 3. Calibrated squared visibilities from VINCI and best-fit 
LD disk model (solid) for GJ845 vs. spatial frequencies. 1-cr 
(dash) and 3-cr (dot) uncertainties are also indicated. 
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Fig. 4. Calibrated squared visibilities from AMBER (low- 
resolution mode) and best-fit LD disk model (solid) for GJ887 
vs. spatial frequencies. 1-cr (dash) and 3-cr (dot) uncertainties 
are also indicated. Error bars amplitudes include both correlated 
and non-correlated errors. 



3. Discussion 

We focussed our study on the low and very low end of the main 
sequence with spectral types ranging from K0.5 to M5.5. The 
stars composing our sample cover a wide range of masses - from 
0.12 to 0.8 Mq - which results in different physical conditions 
affecting their internal structure, the heat transport, as well as 
their evolution. Their atmosphere chemistry is also strongly af- 
fected with the disappearance of true continuum to the benefit 
of a complex and a high gravity stellar atmosphere with strong 
molecular absorptions bands. Metallicity and activity also play a 
role, although we expect it to be a second order effect, at least in 
the near-infrared. In the following sections, we discuss the impli- 
cations and constraints brought by our measurements on stellar 
physics modelling. 
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Fig. 5. Calibrated squared visibilities from AMBER (low- 
resolution mode) and best-fit LD disk model (solid) for GJ166A 
vs. spatial frequencies. 1-cr (dash) and 3-cr (dot) uncertainties 
are also indicated. Error bars amplitudes include both correlated 
and non-correlated errors. 



3.1. Luminosity-radius relationship 



We have first chosen to compare our results to Baraffe et al. 
( 1998| l models, in a luminosity-radius diagram which corre- 
sponds to the observables. Indeed, it has the advantage of avoid- 
ing the inclusion of the mass-luminosity (hereafter ML) rela- 
tionship, for which reliability regarding K dwarfs has not been 
demonstrated yet. Figure |6] shows our VLTI results. Different 
sets of models are overplotted, ie. for an age of 5 Gyr, featuring 
different mixing lengths ( jBohm-Vitense||19 58), L„„^ expressed 
in pressure scale height Hp, that allow to assess the convective 
efficiency, as well as two distinct metallicities : {MjH} - and 
W IH] ^ -0.5. 

[Baraffe et al.| ( |1998| models are in excellent agreement with 
our observations in the very low-mass part of the luminosity- 
radius diagram. The radius determined for Proxima (GJ551) is 
perfectly reproduced by theory. In this part of the relation, stars 
are fully convective which greatly simplifies the modelling of 
heat transport and therefore, our result validates the equation of 
state used by Bar affe et al.| ( |1998| ). GJ 887 is an early type M 
dwarf, located slightly above the boundary of this class of ob- 
jects. The radius determined for GJ 887 is also in perfect agree- 
ment with model predictions. Other measurements from the lit- 
erature confirms that stellar interior physics for this mass range 
are well mastered. 

At the time of writing this paper, there are relatively few radii 
measurements that would allow a discussion in the 0.5 - 0.75 Mq 
region. 61 Cyg A and B radii have been recently determined by 



Kervella et al.|(2008j) and are also well reproduced by theory. 
We note that [Berger et al.| ( |2006| l published 6 radii measured 
with the CHARA array in this part of the diagram. The authors 
claim discrepancies with models at the 2 to 3 cr level. Such large 
departures from theory have not been observed by other studies. 
One may note, however, that 5 of the 6 stars measured by Berger] 
l et al.| ( |2006l l have inflated radii. Those stars were measured with 
the instrument "CHARA-Classic", a recombiner that does not 
include a single-mode filtering. Such measurements are prone to 
systematic calibration errors and indeed, the one star (GJ 15 A) 
which they measured with "CHARA-FLUOR" (instrument with 
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single mode filtering) is in excellent agreement with the models. 
Although a possible explanation, we note that such instrumental 
effect is expected to result in a uniform dispersion. We decided, 
however, not to include those results in this discussion. 

GJ 205 radius, as measured with VINCI, is about 15% above 
models. Radial-velocity measurements on this object have not 
revealed any massive (heavier than a Saturn-mass) companion 
(X. Bonfils, priv. comm.) that would have induced a lower inter- 
ferometric visibility, thus a larger radius. Moreover, GJ205 has 
not been reported to show significant activity (Lopez-Morales 
|2007] l. Nevertheless, this star is more luminous than other known 
objects belonging to the same spectral class and is probably in- 
flated. 

In the upper part of the luminosity-radius relationship, mod- 
els reproduces the observations provided that larger mixing 
length are used (such as L„i^ = l.5Hp and L,„/^ - \.9Hp). 
This part of the relationship may be used to calibrate L,„/^ pro- 
vided accurate observational radii and magnitude determination 
are available. 

Figure|7]displays a zoom on this area of the relation. 

Unfortunately, GJ 663 A does not appear in the tables nor in 
the graphs because of the lack of K magnitude measurements. 
Its measured radius is however shown in table |5] for complete- 
ness. It should be noted that some efforts are needed to obtain 
accurate near-infrared photometry of nearby K dwarfs to tighten 
the mass-luminosity relation and therefore better constrain the- 
oretical models in the upper part of the luminosity-radius rela- 
tionship. 

3.2. Mass-radius relationship 

The translation of our direct measurements into a mass-radius 
diagram requires the use of an empirical ML relationship. We 
used the relation determined by Delfosse et al. (2000) to com- 
pute masses for GJ887 and GJ551, and a recent one, by Xia 
2008|l for the 0.7 to 1. Mq range. This latter study is 



et al. 



based on [Henry & McCarthy] ([1993 ). The ML relationship for 
stars below 0.6 Mq is built on a large number of accurate masses 
and luminosities ( Segransan et al.|2000j ) and exhibits a very low 
dispersion in near-infrared. In this part of the mass-luminosity 
relationship, models reproduce the observations fairly well in- 
dicating that both atmosphere and interior physics of very low 
mass stars are well mastered. The empirical ML relationship 
above 0.6 Mq shows a much larger dispersion than for M dwarfs 
which is related to the modest accuracy of the masses in this 
mass range. The derivation of masses from absolute magnitude 
is therefore less accurate than for the lower part of the relation 
and we adopted an arbitrary error of 5% on masses determina- 
tion. Those results appear in Table |5] Figure [8] shows our results 
in a mass-radius (MR) diagram with results from other studies. 
Five Gyr model isochrones from |Baraffe et al.] (!1998) are repre- 
sented for different mixing lengths and stellar metallicity. 



3.3. Stellar properties 
3.3.1. Effective temperature 

Angular diameter measurements associated with accurate par- 
allax and Johnson photometry allow to derive the star effective 
temperature. UV to near IR photometry have been obtained from 
the literature, mostly from [Morel & Magnen at] ( 1978| l. Visual 
and near-infrared photometry appear in table B.2 We derived 
effective temperatures, by inverting the surface-brightness em- 
pirical relations calibrated by ( Kervella et al.|2004b l. Te/f val- 
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Fig. 6. Luminosity-Radius relationship - Single stars radii vs. ab- 
solute K magnitudes superimposed on 5 Gyr isochrones theo- 
retical models ( [Baraffe et al.||1998| ). Our results from AMBER 
and VINCI are shown as filled circles. We have also included 
other stellar radii from the litterature determined by interferom- 
etry :|Berger et al |(|2006l) for GJ 15 A,|Boyajian et aL](|2008]l, |dil 



Folco et al.|(|2"007K|kervella et al.| ( |2008| l, |Segransan et al.| ( |2603| l 



and |Lane et al. (2001 1 as empty circles. Only radii measure- 
ments better than 10% are displayed. Different models for 5 Gyr 
isochrones are also shown : solar metallicity with L„,u = 1 .QHp 
(solid), L„„ I - 1 .5Hp (dash) and L„i^ - 1 .9Hp (dashdot) as well 
as a metal deficient, [M/H]=-0.5 model with L„ix - l-OHp (dot). 
GJ 205, GJ 887 and GJ 551 that appear in the discussion are la- 
beled. 
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Fig. 7. Luminosity-Radius relationship - same as fig.|6] zoomed 
on the upper part of the diagram. 



ues determined through this method are in good agreement with 
Tgff determined by spectroscopy, both appear in table [5] 

3.3.2. Metallicity 



Metallicity effects have been mentioned by Berger et al. ( 2006] ! 
to explain the difference between four of his measurements and 
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Fig. 8. Mass-radius relationship - masses and radii superimposed on 5 Gyr isochrones theoretical models (Baraffe et al. 1998 1. Our 
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results appear as filled circles. Other long baseUne interferometry measurements come from PTI ( Laneetal. 

et al.|2003| ) and CHARA-FLUOR: Boyajian et al. (2QQ8]), |drFolco et al.| ( |2007l l, |Kervella e t al. (2008) and |Berger et al. '(2006'» for 
GJ15A, all as empty circles. Solar metallicity with L,,,,, - l.OHp (solid), L„„t - l-5i/p (dash) and L„,i^ - \.9Hp (dashdot) are 
shown as well as a metal deficient, [M/H]=-0.5 model with L,„n = l.OHp (dot). Only radii measurements better than 10% are 
displayed. Solar neighboorhood eclipsing binary measurements are represented as empty crosses while OGLE-T transiting binaries 
are represented in filled crosses. Only residuals from long-baseline interferometry results are displayed. 



solar metallicity models. Indeed, the authors claim that miss- 
ing opacity sources in the models, such as TiO, would ex- 
plain the models underestimation of stellar radii for some M 
dwarfs. |L6pez-M orales| (|2007) recently studied the correlation 
between magnetic activity, metallicity and low-mass stars radii. 
Based on [Berger et al.| ( |2006) l measurements, she reach the 
same conclusion. However, no other instrument (PTI, VINCI, 
CHARA-FLUOR or AMBER) could confirm this hypothesis ex- 
cept for GJ 205 (VINCI) as explained in S ect. 3.1. Without th e 
"CHARA-Classic" measurements made by |Berger et al.| ( [2006| l, 
the metallicity-radius diagram for single stars (Fig.|9]l no longer 
shows such correlation. 



3.3.3. Activity 



Close-in eclipsing binaries (EB), for which accurate masses and 
radii have been measured by several authors such as Torres &| 
Ribas| ( p002|l, show sig nificant discrepancies with stellar mod- 
els (see e.g.| Ribas et al.| ( |2008) ). Recently, Chabrieret al.| ( |2007| l 
explained those discrepancies by invoking the reduced convec- 
tive efficiency and stars pots coverage of eclipsing binaries. The 
difference is only observed in mass-radius diagrams (and not in 
the luminosity-radius diagram) because the slightly lower effec- 
tive temperature of EB is compensated by a larger radius, only 
slightly changing the luminosity. This explanation is only mean- 
ingful for EB with periods of a few days implying heavy tidal 
effects, orbital synchronization, and therefore an enhanced ac- 
tivity. This trend is shown on Fig. |8] where EB are represented 
as empty crosses. 
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Fig. 9. Fractional deviation of single stars radii derived by in- 
terferometry from a 5-Gyr, L„„j = l.QHp model (Baraffe et al. 
|19981 l, vs. stellar metallicity. 



However, the same arguments cannot be used in the case of 
single stars. Rotational velocity is an excellent hint of stellar ac- 
tivity. While low-mass and very low-mass EB are routinely char- 
acterised by vsin/ between 7.11 km/s (CU Cnc B, Ribas (2003*) 
and 12 9.5 km/s (OGLE BW3 V38A, [Maceroni & Montalban, 
2004 | l, single stars rarely exceed 3 km/s. We assessed stellar ac- 



tivity for our targets thanks to CORALIE spectra and it indeed 
appears that all of them do not show an activity comparable with 
EB ones. Corresponding log R'^j^ and vsin/ appear in Table 15] 
Furthermore, activity cannot explain radii discrepancies reported 
by |Berger et al.| ( |2006| l since none of those single stars belonging 
to their sample show high activity levels ( |L6pez-Morales|2007[ l. 

Activity cannot be claimed as the source of deviation in 
the upper part of the mass-radius relationship for single stars. 
However, radii of single inactive M dwarfs measured by inter- 



ferometry are in excellent agreement with models from Baraffe 



|et al. ( 1998 1. Thus, discrepancies pointed out by Torres & Ribas 
([2002 1 and Ribas ( 2003| l only concern fast rotating stars, con- 



firming the fact that rotation strongly affects the internal struc- 
ture of those objects. 

4. Conclusion 

Those new results obtained at the VLTI with its near-infrared 
instruments, VINCI and AMBER, allow to better constrain the 
mass-radius relationship for low and very low mass stars. We 
have shown that AMBER is now able of achieving high quality 
absolute visibility measurements provided that at least 3 calibra- 
tor stars are observed and that a careful data reduction and anal- 
ysis is conducted. Using VINCI as a benchmark, those results 
are also shown to be reliable, even if the proposed approach re- 
quires a time consuming observing strategy. Assessment of po- 
tential systematics as well as realistic error bars estimates have 
been crucial to compare our results with models in a meaningful 
way. 

Models are in good agreement with the observations, con- 
firming a correct understanding of the underlying physics of low 
and very low mass stars. The very low-mass regime is almost 



adiabatic and thus constraints the equation of state. A small mix- 
ing length of L„„- 1 - 1 .QHp leads to a progressive underestima- 
tion of radii for early K dwarfs. As expected, the lower convec- 
tion efficiency in K dwarfs require a significantly greater mixing 
length to reproduce observed radii of low mass stars. 
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Table 1. Observation log of stars published in this work. AMBER baseline configuration was always AO-KO-GO (128-90-90)m 
during our observing runs with this instrument. Some targets were observed during the shared risk observing period (march 2002 - 
december 2002) of the ESO Very Large Telescope Interferometer (VLTI). 



Star 


Spectral 


niK 


Instrument 


Date 


Baseline 


DIMM Seeing 


Mean Tq 




Type 










["] 


[ms] 


GJ663 A 


KOV 


_ 


VINCI 


13-05-2003 


B3-M0 (140m.) 


_ 


_ 










26-05-2003 


B3-M0 (140m.) 


0.7 


5 










27-05-2003 


B3-M0 (140m.) 


0.5 


7 










28-05-2003 


B3-M0 (140m.) 


0.8 


6 


GJ166A 


K0.5V 


2.41 


AMBER 


27-10-2007 


AO-KO-Gl 


1.2 


3 


GJ570A 


K4V 


3.06 


VINCI 


07-04-2003 


B3-M0 (140m.) 


0.6 


4 










08-04-2003 


B3-M0 (140m.) 


0.5 


6 










12-04-2003 


B3-M0 (140m.) 


0.7 


6 










15-04-2003 


B3-M0 (140m.) 


0.5 


9 










16-04-2003 


B3-M0 (140m.) 


0.7 


5 










21-04-2003 


B3-M0 (140m.) 


0.7 


4 










27-04-2003 


B3-M0 (140m.) 


0.5 


3 










07-05-2003 


B3-M0 (140m.) 


1.2 


2 










08-05-2003 


B3-M0 (140m.) 


1.2 


3 










09-05-2003 


B3-M0 (140m.) 


0.6 


7 


GJ845A 


K5V 


2.18 


VINCI 


15-09-2002 


EO-Gl (66m.) 


1.1 


2 










16-09-2002 


EO-Gl (66m.) 


1.0 


2 










17-09-2002 


EO-Gl (66m.) 


0.8 


2 










10-10-2002 


B3-M0 (140m.) 


0.9 


2 










12-10-2002 


B3-M0 (140m.) 


1.1 


4 










16-10-2002 


B3-M0 (140m.) 


0.8 


7 










17-10-2002 


B3-M0 (140m.) 


1.2 


5 










19-10-2002 


B3-M0 (140m.) 


0.6 


4 










22-10-2002 


B3-M0 (140m.) 














26-10-2002 


B3-M0 (140m.) 


0.9 


3 


GJ879 


K5Vp 


3.81 


AMBER 


03-10-2008 


AO-KO-Gl 


1.4 


1 


GJ887 


M0.5V 


3.36 


AMBER 


27-10-2007 


AO-KO-Gl 


1.2 


3 


GJ551 


M5.5V 


4.38 


AMBER 


27-02-2008 


AO-KO-Gl 


1.0 


3 



Table 2. List of calibrator stars used during VINCI runs. Angular diameters come from 



Borde et al. 



(2002 1 and 



Merand et al. 



(2004:. 



Calibrator 


Target 


Ang. dist. 


Spectral 


riiK 


0(yo K band 






degrees 


Type 




[mas] 


HR8685 


GI845 


78.3, 19.3 


MOIII 


1.98 


2.01+0.02 


(5Phe 


GI845 


31.5 


G9III 


1.63 


2.19±0.02 


HR8898 


GI845 


12.7 


MOIII 


1.81 


2.31±0.03 


HD 130157 


GI570A 


2.4 


K5III 


2.10 


2.04±0.02 


^Sco 


GI570A, GJ663A 


20.6, 20.7 


K3III 


2.09 


2.04+0.02 



Table 3. List of calibrator stars used during AMBER runs. Angular diameters come from Merand et al. (2004 1 



Calibrator 


Target 


Ang. dist. 


Spectral 




K band 






degrees 


Type 




[mas] 


HD 25700 


GI166A 


9.3 


K3III 


3.16 


1.04+0.01 


HD 27508 


GI166A 


9.9 


K5III 


3.60 


0.98+0.01 


HD 127897 


GI551 


10.2 


K4III 


3.82 


0.91+0.01 


HD 128713 


GI551 


6.4 


K0.5II 


3.49 


0.86+0.01 


HD 130227 


GI551 


6.5 


Kim 


3.59 


0.92+0.01 


HD 136289 


GI551 


10.1 


K3III 


3.54 


0.94+0.01 


HD 204609 


GI 879 


19.9 


K7III 


3.20 


1.14+0.02 


HD 205096 


GI879 


19.8 


Kim 


3.82 


0.82+0.01 


HD 215627 


GI879, GI887 


10.3, 6.9 


K2III 


3.91 


0.83+0.01 


HD 221370 


GI887 


7.8 


K2III 


3.61 


0.90+0.01 


HD 223428 


GI879 


12.4 


Kim 


3.21 


1.07+0.01 
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Table 4. Stellar properties. Masses for M-dwarfs are derived from 



Delfosse et al. 



( 2000 1 while empirical relation from 



(2008 1 is used for K-dwarfs. Metallicity values provided are all directly determined by spectroscopy. 



Xia et al. 



star 


Spect. 


A/I 


left 


Kei. 


Mass 


/TUT 


Kei. 


Derived Tf// 


vsin! 


log R'nf: 


Instr. 


Radius 




Type 




K 




[Mo] 






K 


km/s 






GJ663A 


KOV 










-0.20 


(3) 


4843±134 






VINCI 


0.8I7±0.016 


GJ166A 


K0.5V 


3.90±0.02 


5201 


(3) 


0.877±0.044 


-0.25 


(3) 


5269±35 


0.78 


-4.87 


AMBER 


0.770±0.021 


GJ570A 


K4V 


4.20±0.03 


4758 


(3) 


0.802±0.040 


0.06 


(3) 


4597±101 


1.50 


-4.48 


VINCI 


0.739±0.019 


GJ845 


K5V 


4.38±0.03 


4630 


(3) 


0.762±0.038 


-0.06 


(3) 


4568±59 


1.46 


-4.56 


VINCI 


0.732±0.006 


GJ879 


K5Vp 


4.54±0.08 


4574 


(3) 


0.725±0.036 


0.02 


(3) 


4711 ±134 


2.93 


-4.27 


AMBER 


0.629±0.051 


GJ887 


M0.5V 


5.78±0.03 


3626 


(1) 


0.503±0.025 


-0.22 


(4) 


3797 ±45 






AMBER 


0.459±0.011 


GJ551 


M5.5V 


8.80±0.04 


3042 


(1) 


0.123±0.006 


0.19 


(2) 


3098±56 






AMBER 


0.141±0.007 



- References: 

- For Teff and [Fe/H]: ( 1 ) |Segransan"etaLl p003) ; (2) |Edvardsson et al.| ( fT993) , (3) CORALIE and (4) |Woolf & Wallerstein|p005l l. 

- For vsim and log ie;^^: CORALIE 



Table 5. Derived uniform disk and limb-darkened diameters, and stellar radii. Parallaxes are from Hipparcos ( van Leeuwen 2007 1 



Target 


Instrument 


parallax 






Limb darkening coeff. 


duD 


Sld 






[mas] 








Kband 




[mas] 


[mas] 


GJ663 A 


VINCI 


I68.54±0.54 






[0.79, 


-0.56, 0.43, 


-0.15] 


I.253±0.025 


1.282±0.026 


GJ166A 


AMBER 


200.62±0.23 


2.39: 


t0.02 (a) 


[0.81, 


-0.53,0.39, 


-0.13] 


I.405±0.038 


1.437±0.039 


GJ570A 


VINCI 


I71.22±0.94 


3.15: 


t0.02 (a) 


[0.86, 


-0.52, 0.37, 


-0.12] 


I.I47±0.029 


I.I77±0.030 


GJ845 


VINCI 


276.06±0.28 


2.18±0.02 (b) 


[0.86, 


-0.53,0.38, 


-0.13] 


I.834±0.016 


I.88I±0.017 


GJ879 


AMBER 


I31.42±0.62 


3.95: 


t0.08 (a) 


[0.86, 


-0.53,0.38, 


-0.13] 


0.750±0.066 


0.769±0.067 


GJ887 


AMBER 


305.26±0.70 


3.36: 


t0.02 (a) 


[1.61, 


-2.35, 2.00, 


-0.68] 


I.284±0.031 


1.304±0.032 


GJ551 


AMBER 


771.64±2.60 


4.38: 


t0.03 (c) 


[1.94, 


-2.80, 2.39, 


-0.81] 


0.990±0.050 


1.01I±0.052 



- References: (a) |Morel & Magnenat|(T978l l; (b) |Mould & Hyland|(T976l l; (c) 2MASS l |Cutri et al.|20"03l l and (d) |Segransan etd^jpOOSl l. 
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Fig. A.l. Calibrated squared visibilities from VINCI and best-fit LD disk model (solid) for GJ570A vs. spatial frequencies, l-cr 
(dash) and 3-cr (dot) uncertainties are also indicated. 
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Fig.A.2. Calibrated squared visibilities from AMBER (low-resolution mode) and best-fit LD disk model (solid) for GJ551 vs. 
spatial frequencies, l-cr (dash) and 3-cr (dot) uncertainties are also indicated. Error bars amplitudes include both correlated and 
non-correlated errors. 
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Fig.A.3. Calibrated squared visibilities from AMBER (low -resolution mode) and best-fit LD disk model (solid) for GJ879 vs. 
spatial frequencies. 1-cr (dash) and 3-cr (dot) uncertainties are also indicated. Error bars amplitudes include both correlated and 
non-correlated errors. 



Table B.l. Online material - Angular diameters and linear radii of single low and very-low mass stars from the literature. They are 
respectively expressed in mas and sun radii. 



Star 


Spect. 


[Fe/H] 


Ref.[/re/H] 


Instr. 








Ref.,„,j,.y 


crDra 


KOV 


-0.21 


(a) 


CHARA 


1.224: 


t 0.011 


1.254: 


t 0.012 


0.778±0.008 


(4) 


HR511 


KOV 


0.02 


(a) 


CHARA 


0.747 : 


t 0.021 


0.763 : 


t 0.021 


0.819±0.024 


(4) 


eEri 


K2V 


-0.06 


(a) 


VINCI 


2.093 : 


t 0.029 


2.148: 


t 0.029 


0.735±0.005 


(1) 


GJ105A 


K3V 


-0.07 


(a) 


PTI 


0.914 


±0.07 


0.936 


±0.07 


0.708±0.050 


(2) 


61CygA 


K5V 


-0.27 


(a) 


CHARA 






1.775 : 


t 0.013 


0.665±0.005 


(5) 


GJ380 


K7V 


-0.03 


(a) 


PTI 


1.268 


±0.04 


1.155 


±0.04 


0.605±0.020 


(2) 


61CygB 


K7V 


-0.39 


(a) 


CHARA 






1.581 : 


t 0.022 


0.595±0.008 


(5) 


GJ191 


MIV 


-0.50 


(b) 


VINCI 


0.681 


±0.06 


0.692 


±0.06 


0.291 ±0.025 


(3) 


GJ887 


M0.5V 


-0.22 


(c) 


VINCI 


1.366 


±0.04 


1.388 


±0.04 


0.491±0.014 


(3) 


GJ205 


M1.5V 


0.21 


(c) 


VINCI 


1.124 


±0.11 


1.149 


±0.11 


0.702±0.063 


(3) 


GJ15A 


M2V 


-0.45 


(d) 


PTI 


0.976 : 


t 0.016 


0.988 : 


t 0.016 


0.379±0.006 


(6) 


GJ411 


M1.5V 


-0.31 


(a) 


PTI 


1.413 


±0.03 


1.436 


±0.03 


0.393±0.008 


(2) 


GJ699 


M4Ve 


-0.15 


(a) 


PTI 


0.987 


±0.04 


1.004 


±0.04 


0.196±0.008 


(2) 


GJ551 


M5.5V 


0.19 


(f) 


VINCI 


1.023 


±0.08 


1.044 


±0.08 


0.145±0.011 


(3) 



- Ref. for [Fe/H]: (a) |Soubiran etaTlpOOSl i; (b) |Mould & Hyland| ( (T976) ; (c) |Woolf & Wallerstein| ( |2005) ; (d) |Bonfils et aTlpOOS) !; (e) |Edvardssonl 
[etal. (1993). 

- Re f for interferometric m eas urements: (I)|di Folco et al.|p007| l; (2) |Lane et aTlpOOl^ ; (3) |Segransan et al.|p003) ; (4) |Boyajian et al.| ( |2008l l; 
(5) |Kervella et"ar]p008l >; (6) |Berger et aL] ( |2006| l. 

Table B.2. Apparent magnitudes of stars included in this study. The uncertainty adopted for each apparent magnitude value is given 
in superscript. 



Star 




















GJ166A 


5.69»''2 


5.25»"2 


4430.02 


3 72""'" 


3270.01 


2 9^0.03 


2460.01 


2390.02 




GJ663A 


5.670-02 


5.18'"'2 


4 320.02 


3^62002 


3.18''«2 










GJ570A 


7.88''»2 


6.82°''2 


5Jl''-"2 


4 720-02 


4 180.02 


3.82°-02 


3 27O.02 


3.15«-''2 


3.1l''-''2 


GJ845 


g 740.02 


5.75»"- 


4.690.02 


3'81''»2 


3.250-02 


2.83"-02 




2 lg0.02 


2 12002 


GJ879 


8.61''»2 


7590.02 


6.490.02 


5.540-02 


4950.02 


4.53°-0* 


4.00'"" 


3.95«'"* 




GJ887 


9 990.02 


8.83'"'- 


7.35''-"2 






4 200.02 


3.60'"'2 


3.36«-''2 


3.20''-''2 


GJ551 




13. 02" "2 


11.05''»2 


8.68''''2 


6420.02 


5.33»-«2 


4730.02 


4.36'"'3 


4040.02 
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